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PREFACE 


If experimental work is to be accomplished involving 
the use of a nuclear reactor, accurate knowledvemecween aad 
bhEe density of neutrons at various locations withing ene sce. 
actor is a necessity, whether the reactor be crijmica teen 
Sue ri FLical. 

The néutron density distribupion in aeons oe on 
can be estimated mathematically by the use of SseVeroisve ch. 
Higues including transpert theory, multigroup dittuo tema 
Ory, and single. eroup dismusions theory. Of these various 
Mathematical theories, the single group diftivscitommcea. 

Omly one amenable to@hand caleulatiens. Untomtunavely shew 
ever, diffusion theory is alse the Weast accwrat wees 
three named, althougheit does provide fageecode tine appre. 
MmepetOn to the actual neutron density distribution, wives 
Eeacuor. 

In order to accurately determine the neutron den- 
sity distribution within a reactor, experimental measure- 
ments must be condueted using one of the various methem. 
evollable vror ss UcimemeasUt ec meme cr One of the most "ceo vem en 


methods is that of foil activation measuremenitce. 





Control rods made of materials which can absorb 
great numbers of neutrems are used = to emtrol them@meutron 
@istribution within regetors. ~in critical Syetewoee oie 
Gontrol rods are, of course, necessary for ‘SsSatemeand paar 
Operation of the reactor, 

The Universitg,y of Texas mMamnutains a suberieweca 
reactor which is used for teaching purposes and for research. 
whe purpose of this thesis is to aceurately determi. 
effect of a pure cadmium control rod on the neutron deme, 
fis ttt bUuGLOn in thse reac cor. 

The neutron density distri~ution is deternmicdy.., 
fOoLtl activation measurements both With and without theweacd= 
meme cOonbrol rod. Single group diffusion theory is used to 
predict the distribution mathematically, and the theezve Gia 
distributions are compared with the experimentally determined 
ones. 

The experimentally determined neutron density sai 
bri pution without the control rod is compared wich neo 
tribution with the cadmium control rod in place Wier 7 aa 
elearly show the effect of the rod on the neutron density 

The author wishes to express his appreciation to 
Dr. R. N. Little for his advice end guidance dunia 


Supervision of this research and to Dr. J. AP vscentanewan 





his many helpful suggestions concerning the writing of the 

thesis. Sincere gratatude is also due Mee or Ow Maethur “and 

fie, C. J. Opahn for’ their patient, expert iiyande eres comm 

summing assistance with the equipment used in this experiment. 
Since the writing of this thesis marks the end of 

a long, hard, three year struggle, the author fee eea tae 

snould at least Nave a dedication, Therefore tii Sow ieee 

is dedicated to the person who has provided the €ncoUreagenenm: 

and the understanding which is so necessary in this@i yp emer 
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CahiA PP TER it 


THEORY 


A. One Velocvty Croupebiitusion Theory Approxima mem 
Rods 


WieWOout Control 


In the single velocity group dititieamon cheer sue 
assumed that all neutrons in the system being studied have 
ene same velocity, which in the present case 15 2200 meu. 
mer second, ene velocicy ef thermaleneurtrom., iit eG ie 
neutron population is described in terms of the neutron 
density, n(r,t), where 

n(r,t) = number of neutrons of speed v per unit 

volume around space point r at time Jt, 
then the neutron flux, ¢g(r,t), is expressed in terms of 


men.) by the usual definition 
g(r,t) = vn(r,t) 


maere, in this case, all valves are for wherna ememmom= 
and v is constant and uniform throughout the reacvorn. 
A verbal statement of the neutron balance jequawwen 


hor @ differential volume element, dV, inclu@es toe aie 


Pi V. Méeghreplian and D.” Ko dolmes, Reaewomm ace acess 


(New York: McGraw-Hill Book Co., 1960), Chap. 5. 


alt 





@ 
account for the net gain or loss through the boundaries of 


ecu y OLUmMe wellomen > seams 


Or =the meutron 
density in dv 


fa OL citron 
about r 


duced in dv sorbed in dV ing fromeeasy, 


HhuUMnber ot os ous | Re number on | 
DeEeuni ts time 


Pose pro- neutrons ab- neutrons e€s¢ap 


per unit time per unit time 


Mas statement, in mathematical form, becomes thevequatuen 
ee Vr iy Gd J emi gneve erica 


where S(r,t) is called the source function and denotes the 
Momoer Of nNeutroius produced per Whit time Andy eee 
term &£(r,t) is called the leakage function and denotes the 
itt LOSS Of Neutrons from dV duce CO “ditiucion hwo a. 
boundaries. 


in diffusion theory the leakage term becomes 
OG Gaas), a a DV“ ¢ Ee) 


and the neutron balance equation becomes, UnderesG-a0 esa 


Conditions for which nine hb) iseconauami., 





~ DV“g(r) + zo (r) = S(r). 


in this equation D isvealled thevdiffusion coer teeren pean 


is defined to be 


cer ile 
S on pe 





Z (cos 6). ] 


13(@.a0rs mb veg 


where A+) is the transport mean free path for | thes cue 
eontaining the volume element dV for which the edtjaticn aa. 


aire Weoly Be LS the macroscopic absorption cross Secrweae.e. 


a 


this medium, and S(r) is the source term. As indicated above 


it 


rv Sm cs OE EE ee ee 
OIC 
Le cvamn ats BECeos IE ae | 


where — LS the macroscopic Scattering Crosseeecomonmor ta 


medium, (cos @) = 2/3A and is the average scattering angle 


avg 
for neutrons involved in seattering collisions Vr cee ee onn 
the medium. 1D IS el avowesthit lei eke fell 2, 

total S a 


The net neutron current through a given area lo- 


eated at r is given by the vector equation 


I(r) = -DVg(r). 





4 
tie partial neutron euirrents through an aieewear oc dciinced by 


the equations 


j,(r) = = ie i 2 \Vg(r)]|cos Qo 
and 
j_(r) = = g(r) + : lvg(r)|cos 6, 


where oe is the angle between the unit vector normal to the 
area and the unit vector parallel to the vector Vg(r). The 
@ifference of the partial currents is equal to) (een 


omer words. 
I(r) = 3 .(r) - a_(r). 


The boumdary conditions which must be used in obtaining a 
Soluweronw ror thewabove diffusaon theeryieutrenm Datemce sean. 
fom are as follows: 

(1) g(r) must be finite and nonnegative in a1] eevons 
of the medium in which the diffusion equation ap= 
plies. 

(2) At any interface between two media witheatt tone 
diffusion coefficients, the partial eur Vem sees 


bie eon cine s = Giolla eer Giere st oem Lipebmere Ts 





5 
no external source at an interface ecient tic Wewiron 
flux, g(r), and the net current, J(r), must be con- 
tinuous 24, tthe an ter tace. 

(3) The flux g(r) must vanish at some finite distance 
outside the boundary or else the partial current 
from the area outside the boundary must be zero. 

In diiiusion theory, the distance be yo mdg Gg emaes 
actor boundary at which ¢(r) = 0 is called the ex- 


trapolation distance, d (= 0.7104 Lee) |< 


Two assumptions were made in this treatment con- 
eoonine the reactor vassembly without che contro meaeas iat 
Was assumed that the air gap along the central axis of the 
eore did not exist. PRis air Gap actually represen rcameus 
approximately 1.2cO per Cent of the vVetume of tne tome ee 
Giks ASSUMPLCION was Not considered unreasonable. It was 
also assumed that a point source of neutrons was located at 
ime center of the reactor, This point source is treated as 
a boundary condition in the following mathematical develop-— 
men tc. 

In the case of finite cylindrical 2seome cue, eco 


frtiusedon equation can be written 


7G DV“3(0,2) se Tazo, Z ) = Son 4), 





ae 
where S(p,z) = vEpp(o,2) for a multiplying medium and is 
Zero fOr 8 NONMULDI DE yams medium, 2 is the macroscopic 
fission cross section for the medium and v is the average 
number of neutrons resulting from a single fission event. 


Li both a@@eone wand a reflector are imvie ew aye 


diffusion equation must be written for both restons sete oyvs- 
DV fel0,2) + vEe Plo.) = Ea Php, 2) = Q 
DaV*BR(0,2) - Za alo,z) = 0 


The solutions to these simultaneous equations must satisfy 
all the boundary conditions previously stated. 


The aboves@equditaeusmeanmalso be written 


Vp (on ) eaeBe acne 


and V-brp(p,z) - KRP (p>) = 0 


e = = 2 = Cio 
where Bo SNe EVD, and Kp = Zap! PR’ In 1eUi Gio Waiceinns 
fusion theory a term called diffusion length, Ey eee dewrted 


such that 


D 
2e- 2 
Go = ¥ 





Using this definition eseim the terms a00yemwerve> 





Ke —= pe 

R fe 

R 

Vea k - li 
f 

and Be oS ee gett he 
c Pg oY Le 
CG ac Gs 


oeGe kK. = 


vie /Za_ is the infinite multiplication constant 
Gs Cc 


for a medium composed of the same material as the reactor 


wore. 


In cy lint@rircal coordinate the vo Opera COT =a 





so the above diffusion equations become 








19 | = 
a ac ae se(0,)] 2 Y Sa £.(0,2) An Ia) 275) 9/4 | 
193 | o- 

and a 2 an dao.) A ae pp(o,z) Kp plo,z) 


Doth of which can be Solved by the 
variables. 


For the core it is assumed that 


Fale Ufo 249, = Re Govze (zy) 


| 
© 


| 
) 


method of separation of 





8 


and the resultant general solutions towthe separated ditf- 


fucton Gduavion fOr time "core ane. at Be = ke + Ke, 
Yo) = he Sim kee as Cl Cosma 
c Z Z Z Z 
= Ji + k 
yao) See se) Ce) 


where J,(k,0) and Hee) are zero order Bessel functions 
@i the first and second kinds, respectively .| [hee eumcianen 
Yo(k,o), which approaches negative infinity "aso apereachnes 
Phe value of z2éro, becomes more negative for Smaliievailuee 
ohn ao) than the function Jo(K,oe) becomes positive. walhuc 
bumee the boundary Conditions Stipulavesienat b(0,z) can 
Meyer be Negative Cy Must beetakem to bemzZero. Weamice 

san K>zZ becomes negative when Z becomes Nescative tis een 
would violate the physical requirement that 6(p,z) be syn- 


metric about the z= O plane. Thus; A must alco lbeur aia 


Z 


wo: be zero. 


Hence, the solution for the come Trecimomaia 
bya (og) = A oJ g (Kae) COcwK 2 


For the reflector it is also assumeaw vat 





po(0,2) = RA(o0)Z,(z) 


and the resultant general solutions to the separated dif- 


MicltOnmeQUaulOn [Cr taerrel lee Comware., it Ke = Ke = Ke ; 


Vi Stes ae 8 COs 4 


we) 
oe) 
[2) 
lI 


i el 2 + C Ko (Kk) 
where eno) and Ko (t,0) are zero order modified Bessel 
bummer ions of the first and Second kimds respect voine Here 
meain, the condition of Symmetry Of the neutron tite a2 voy 

Z —- 0 requires that A, be taken to be zero. 


Hence, the solution fOr the cea leeiOn sec aie nee 


Pplo,z) = [ApI, (6) a CaKy (Qo) ] cOS KZ 


It po andeq, are defined fo be tie cxirapo lies 


O 


madius and height, then the flux must be zero at Tremp amir. 


where 96 = O59 Or where Zz = Z). This boundary condi (tone: e— 
quires that Pp(o0,z2) = O when 9 = pg Or when Z = Zo. 90 
Ke Zy = t/ 2 since coc =a One cains Kk, = t/ 22. 


This boundary @ondition alse requilreewunas. 





10 


[ApI,(« 0) + CRaKi (Ko) J = O° ~ when aos Oo: 


Thus 


Mie Equations for tee flux in the reactor thus piecemeal 


summary , 


Pia ye) = AgJg(k,0) cos k,z 


Ko (« 09) 1 
ép(o,z) = Cp Ko (0) fa TolEsbo) 1s aren), cos 22. Z 
mere Z, and p,.ame the extrapolated hementeyand radiu2 aoe 


epeculvely. 

At the interface between Ghe core anew tne wre 
mbeeuOr, the boundary conditions require that tie sili sane 
Peemner Current be CONTINUOUS aeCr@ss the anmtecriace. lem ere 
the current i8 a vector quantity, this conditions. cape. 
baat the components of the current vectors De (eCitenve e 


daus, at the interface denoted by height 2, and Vadii a. 


Gee | = Pplo,z;) 


and To (0,21) = TR 0,24) 





ala 


which leads to the equations 


Kein om) i 
AgJo(Kpo) COc wie. 2 on oxo) - TEST Fol toe) 608 Pas Abt 


OD 


and 


DekeAgdo (Kao) sin k,z; 





Ge Stee) 
a Oo O Tt 
re non oP)" ta es a Bay 4 


Dividing the second of these equations by the first gives 


the equation 








k So 


R 1 
tan Ko 24 = ae tan ao 


Z 


mien can be solved by numerical methods ctoMmdevcrminec ek 


and Z. are cOmnputed wor meas 


after the values of Dp, D a 


on, “i 3 


ured for the reactor being studied. 


Once the values for k,~ and k> are KiG@Wi eu co. 


Z 


responding values of Ko and kK, can be computed from the re- 


O 


ieab1Oons 


Be = Is + ke and KR Se ee Re 


ia 


i 





12 
since Be and KE depend on the physical constants of the re- 
He bor and Chese  COnsGants aGeukniown. 

The constant A. is related to the constant Ce by 
tae iInvcerface boundary condition that the flux must oeecon 
tinuous at the interface denoted by radius op, and the height 
wee hus pa(os52) = polo, ,2) or, if the point where z= 0 is 


Ronen, 


av J,(k.9,) 
aS as oo. 
Fo(epes) = er, Fo(noer)} 
Oo ialae, 
Hitch will give Cp as a eonstant multiplied by Ac. 
ive (CO ls ue mt oe can be evaluated by the boundary 
Seemdition of continuity of net neutron current al yene io 
face between the external neutron source and the core. If 
Somparison of theory with experimentally determined ip 
Wie bribubtions is desired, then the maximum value ot rne Vie 
ean be used for Ae 
The values of De, DR, Be, Ky, Kz, KR; Ky, Kos and 
Cp for the reactor used in this experiment are computed in 
Appendix A. 
If the neutron density distribution, n(p,2), 


rather than the flux distribution, §(p,z), is desired, the 


the constants A, and Cp are simply divided by the thermal 





ES 


meutron. velocity (27165 O00 cm/sec) to obtain the equations 


= 1 
n (0,2) = AQT, (0) cos kz 


Ke Uk oie) 
O O 
nplp,z) Ch Ko («,0) = rep tePO 5 (x90) feos Ke 2 


Z 
T(r 6) 


where 1 = A./ Vin eulaiie| Ch = Cp/Vin- 


Pee One Velocity. Groupe iis von theory 
Approximation with Control Rods 


Oince the diffusion theory differential equations 
are the same for the core and reflector regions regardless 
fia tne presence of a control rod, the general score cw mic 
miese CQuations remain the same as for the easecavnewmenew en. 
meol rod is present. The effect of the control rodeisieec— 
Seounved for by the boundary conditions used to Cvaluavc mr] 


eeonstants. 


The general solutions for the diffusion eq ione 


ale 


peo; 2). = [A.JQ(K,0) + Co¥o(k,o)Jeos Kee 
Pplp,z) = [AR tS Re) “ CrK,(« 0) Jcos ee 
a 


R. V. Meghreblian and D. K. Holmes, hegetor 2.) 
sis (New York: McGraw-Hill Book Co., 1960), Chap. 11. 





14 
where Be = Ke + KC and os = eA = KO as before, and Jo (kyo) 
and Yo (k,0) are zero order Bessel functions of the first 
and second kinds, respectively, and I,(,) and Ko (0) 
ame zero order modified Bessel functions of the fies eee 
second kims, respeceilve ly. 

The boundary conditions which were applicable for 
meme Case with no control rod@remainvapplicable, and forge 
eise Of a black control rod@= fully inserted alongs the vaxas 
or a cylindrical reactor, tne additional boundary condi 


that Po = UW al the point Wiereso.— Denis Dewampecca- The 


O 


Malus bo 1S defined py the relation 


mere bis the radius of the control rod and e« is the x 77. 
polation distance into the rod. 

When a black control rod as present, stnewt i eee 
comes zero for 9 >O, and so the term Yo(k yp) need not be 
eliminated from the equation for 6,(p,z). The boundary 


Soudition that ¢.(b,,2) = © yields the relatiensiag 


Agdo(Kybo) = - Ce¥o(kybo). 


ifs relation can be used to eliminate ©) anda ea 





equation 


where A,»k and Kk Saree bau Dendevermined, 


0? 


The boundary conagdi tion thia u dp(o,z) = 0 80 view 


teapoOlated radius, p24, and heen, 2 yields, as fer awa 


Or? 


mee Of NO COntrol rod, the rmewarionshtps 





I 
Ke ~ O¢ 
Keay Gheaoy 
and Ap ies Cr Tie a 
en alee 
where Ky can be determined from the definition KR = ee _ ke 
when Zo Ls weno wi, 
in summary, the equations fOr the 7 luweeomeno wer 
written 


Jo(Kpbo) 


Crone) be [Folge) = ¥, (KB) rot 0)| cos k,z 


Ko (K,9) 
iene) 
oo 


li 


pplo,z) 07 [Ko (ego) Zo 90)| COS Jka 


je) aa = K are known. 


miere kK we 0 


pers 





ead 


dG 
The boundary vcondition for contimue, Of flux sand 
net neutron current are used to determine the remaining con- 
stants. 
At the ingverface denovedsby heleht 2a te ema, 


the relations 


p.(9,2;) = Pplo;2;) 


and Jo (0,23) = Jr (e524) 
must be true. These relations again lead to the equation 
cos KZ, : cOS KZ; 
Dake Sin kez; Mews | Sabie: Weve 


which can be written 


DR if 
ken Sean: Ky Zt = ee tan 


e Day 27 EMVags 








oo 


and can be solved by numerical methods in order to determine 
Ko. After k, is determined, ee can be found freon) vaew 
PO 1B = ke “+ rae which has previously been defined. 

The relationship between Ac and CR is determined 


by the boundary condition that ¢,. = gp at the interface 





ule 
denoted by radius op; and height z. This condition yields 


tie relationship 


Dee) 
nea) aa eee) Mie alos? 


Ko Ea Oo) 


CR = Ag 
K_ (kK . == eee ~ RD 
mi the z= 0 plane is used for the evaluation. 


In a Similar manner, the boundary condition uae 


Cc 
O 


te- UO plane can be used to determine A since this. ei 


J - Jp at thie interface Gdenoved by radius o> sane 
O 


two simultaneous equations containing A, and Cp. In order 
to correlate theory with experiment, A, can alternatively be 


taken to be the flux observed along the midplane at the dis=- 


fomec 9 Lor which 


Jig (kaon 
0 (%90) = 7, (k,bo) toltge)| ae 


The values of all constants except A) aude ae 
mime same as when no control rod was present, The values of 
fm, and Cp for the reactor used in this experinicni aa ear. 


c 


puted in Appendix A. 





As. 
If the meutron density distribu ou, p47 jes 
ected urather Chan teem rh xy di Ci bu trom b(o,z), tie Com 
stants A, and Cp are divided by the thermal neutron velocity 


(220,000 cm/sec) to obtain the equations 


ad Gl elev 


Ba fTolg0) - T. Edo) Yo(k0)] 08 kee 


leery 4) 


ik ce) 


oh] (moe ; Tern) 1o(*p2) £08 Se 


6 eek e570) 


np(o,z) 


where Ai = A./vsy, and Ca = Cp/vip. 
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NEUTRON FLUX MEASUREMENT BY THE ACTIVATION METHOD 


in. General 


The activation method provides one of the most con- 
venient and most accurate techniques for measuring neutron 
16 JL ip ge itis metiod is based ™@on tre fact that many -elememre 
become radioactive when exposed to a neutron flux. ror. sar 
sual elements, this activation occurs “only ToOrvae narrow 
range of neutron energies. By measurement of the induced 
meaivyitcy, the neutron flux responsible for tae acum ome 
can be determined. 

To measure. a flux; a echimwetoil of Une os eer mua ee 
Maverial is first irradiated for a certain time interval at 
menaesired location. The foil is then removed  erom ene 
fmeutron riux and the induced activity is measured “inwan ap 
propriate counting device. The relationship between the 
measured activity and the neutron flux must inelude "correc— 
Giron terms for such things as counter efficiency ena, te. 
eon of neutron flux by the foil being@used (oO perio wiae1e 


measurement. 


ee 
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In order to be suitable for this type of measure- 
Ment, the foil material sh@uibad Nagswe a large neutron abso;p— 
Puen cross section Tor the emerey mange to bewstudaed war ltso., 
the decay scheme of the radioactive isotope produced by the 
Meutron activation should be simple enough for eas, measure— 
Hem@imand interpretation, and the halt? (ite for decay soie cite 
product isotope should be such that irradiation and count- 
mie periods are Mot excessively lone. 

For thermal neutrommdietectron, the mos teireduvenman, 
used Sisters which have strong resonance absorptions and 
Vaetd induced radioactivity of suitable half lives and Taga 


mronscvare in, Rh,eAg: Au land Dy. 


Be Use of Py sprosium ois  Toreiiermor 


Neutron Detection 


The dimensions ofthe suberi1ti cal sreac Omer et 
bis project made necessary the wse of a Van de Gragif ac— 
welerator to produce the external source ef nevtroems meecs— 
Peary for a subcritical reactor (see Chapter Ii] ) meme 
meee operational characteristics=of the acceben geste. ode 


was possible to maintain a constant source strene vito ee, 


Sc. F. Bonilla, Nuclear Engineering (New York: 
MeGraw=Hill” Book Co.4. Ime secon Dp emece e 





ay ah 

SiOmie per LOdSs Of “Tiine and Jenence, the toile vee tote ma— 
terial used had to be such that maximum induced activity 
could be reached in the shortest period of time feasible. 
Of the foil materials (In, Au, and Dy) available, dysprosiun, 
which has a large cross section for thermal neutron absorp- 
tion as well as a low cross section for epithermal neutron 
beer puvlon, best met this requirement. However, dysprosium 
does have one disadvantage in that approximately twenty 
hours must be allowed for adequate decay of the induced ac- 
Pavity prior to reusing a particular foll for another snea.— 
Peeement. 

Preliminary calculations showed that ot scien ca. 
seal naturally occurring Brsotopes. of dysprositumevaleh wood 
meoltlt in radioactive products when irradiaved )) sueutrouce 


164 


my the isotope Dy need be considered. Contr pl cioms 


co runduced activity from Weutron absorption by isovopes 
ieee 2° and Dy+°® are negligible due to low percentarcauor 
abundance and small neutron absorption cross sections as 


well as the long half life of decay of the radiodactive prod-— 


Ley Sg 


mets Dy and Dy The contribution to measurab tow a 


tivity by the transformation of Dy*©®° into (the Trad meeera 


tec 


asotope Dy was determined by calculation to be less than 


gax 1Om> per cent: These calculations were confi emeamp, 





- 


aa 


experimental determination of the half life of decay for 
EHe activity induced in mMatural G@ysprosium byenreulr one ac— 
tivation. The observed half life was exactly that of Dy‘t®°. 

Thus, onlys the activity anduced throueh aie naaie® 


164 need be considered in con- 


absorption by the isotope Dy 
nection with neutron flux measurements. The reactions in- 


volved are: 


Dy -Y 
D 164 165m ae 
y + n >~ Dy rage 
HOt or 8B +yY 
Dye oem eelvege Silica =. ang mepe 
165 


iiwe decay scheme for Dy is shown in Figure sei 
The general differential equation for formation 


Sieraagioactive product nuclei by the neulvon ea eee wel. 


meron 15, in word form, 


Rate of accumula Rate of produc- Rave Di eloss 
Laon “Of radio— =|] tion of radio=— |= omer adie: 
Seuive muclea active nuclei active nuclel 


im the case of Biya a two differential equations are neces- 


Les 


sary since Dy?°* can be transformed into both Dy and 





p2 


Bl 


e5) 


Iie sae (1.25m) 
O.108 Mev 
nA 8 decay 
I7.6% 2.4% 
2 5 Sa) 
ps = 
Mev 


G5 05 — Mev 
7 f 


ee |. 
de O55 
es fe Hip 
. GnS5e 
ALLL oes 
Mece ot || 
0.095 





Hoss 
: 165 
Figure II-l. Decay Scheme for Dy (from 


Tornau*®) 


*Ruidiger Tornau, Zeltscehritt fur Shy coe eee 
OAs ALG ey. 
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ilgeysy ddl 


Dy by the neéutroneapcorovlonrrecac lion. These equations 
ae 

oD. ae) 

pee OM ie pA om z 
(1) aT > oe DaPth ~ 5m25m(T) 
( 2 4-Sp 4 976 X._D__ { 
) aT = ees ae 5m 5m‘) 

- = Sie 
he De (T) ore Do(T)Pi, 


The notation used in the above equations is as follows: 


164 


Da =—"Humber of Dx ig Ole Jie alee Oe ers (Sal ae (assumed to be 


Come want) 


aon s 


ea, 
on 
— 
F 
— 
; 


humbier Of TDy mueled present av tame wl. 


ie Sl] 


Do,(T) = number of Dy Auclei present au enc mle 


AL = TPRaead es Giem Gaeine 
Den =Jthermal neutron fTlux duvimey irradia viom 
Tae = thermal neutron absorption CEOSs, seceu neuen 


Hevive tion cima. amepouD alae 


as = thermal neutron absorption cross ‘seeUleneier 
Z : : 164 165m 

aAGCCiLVablon Oia s tO eb y 
ony = thermal neutron @bscorption Cross Sect tenwenon, 


Aetivation or Uy. mom Dy. 2 = 





cs 


» = decay comstant, ter Dye 


» = decay constant for Dye 


5m 


Preliminary ecaleculations showed that the eilicertsor 
Ee activation of Dy > to Diy is completely merci eaeier 
SO tne solutions to these two simultaneous differential equa-— 


muons are found, by standard methods, to be 


4° SMH pb 

oO 

_ ny 4° th 

DoT) = oe 1 - exp (-A5mT) 


and D.(T) = Bf - exp (-050)] + cent (eee) - exp (Ag) 
4-5 4-5m 
2 oy ae 
where B= x bth 
Ss 
4-5m 
7 Ooo C., D, 
and Cc = ———_—_+__| Dth 
(Asm 7 5) 


im determining these solutions the boundary condition that 
both a) and Det) were zero at time T = 0 was used. 


If the foil is removed from the neutron i a. 


Les 165m 


meme I, the number of Dy and Dy nuclei present at any 


time t thereafter is determined by the solutions to the dif- 


ferential equations 





C6 


GHDeen et) 

om = - Nep~Dey(t) 
cE Oiler) 
abet) 


In this case the boundary conditions to be used are that 
Devt) = De, 2) and Delt) = Del?) at t = 0, where = teutime 
acer removal from Che neutron flux. 

The solutions to these simultaneous ditferentrer 


equations are 


Dee = Dgn() [exp (rsnt)| 


OVeers uae 


and Det). = Pen 7) [exe (- X5,t)- exp (- 5) 


+ Dg (7) [exp (- 5t)} 


DIince the activity at any time UC ls Simply some 
Peoauct of the decay constant and the number of Mucous — 
ent at that time, the following™equations give thewac Umer, 
G@>o a function of the time after rémoval of the Gece eee) 


iom the neutron ilux: 
Aew(t) = AgmDem it) 


Ac(t) = AcDe(t) 





a7 
THUS .07t “thes aetryi ty 1S  COUl Poirot. om. fC eis 
Ponce incerva tl. ty to to» the number Of radwocdeu ce  Uieioe — 


Spacious, N, Oc¢curning during tais time wil eie 


We e, 
N = A, (t) ie A(t) at 
ty ty 
bo We 
= Lar) dt + WD. (t) ae 
ty ty 
where N = total number of disintegrations occurring between 


iemes th and to Evaluation of this integral yields 


-- 


0.976 d, 
im |D.( 2) - (ee = ey Da) ext (- A,t1) - exp (- gta) 


r x 
5 5m 


A he - Aon 


| (CS) eee Dea exp ae) ae crse ae sate) 


this. equation for N appears to be extrenc ain. 
Bersome and unwieldy, but after the constants arege waa 
amd the time intervals to be uséd for irradiation sanaeeeia 
ing are derided upon, the above equation, reducocmrem a 
eonstant X m Pin where m is the weight of therde vec vom a! 


in grams. 





ao 
trial -ruds wand. caleubations Wedge iots wo reve hous 
weights indicated that an irradiation time of five minutes 
(T = 0 to T = 5) and a counting interval of 10 minutes 
(t, — ONO te = 25) would give very accurate results. Evalu- 
ation of the various constants needed for the above equa- 
ions Was carried out and the resubis are slanuia lode 


ioe “t-te Usime these eonsctan espero weauavl ony ter wie. 


Geeces LO 
N = 0.6402 mpiy, 


Mere m= foil mass in grams and bin = neutron flux in 
Mmewtvrons per em“-min. After converting, Whe trl luxe avoe mel 
Meems per Square Ceimeimeter per second, viewedua pilot. 


thermal neutron flux can be written as follows: 


J neutrons \_ N 
th\ em? sec 38.412 m(grams ) 


Pee De cerminationwoL [hermal Neutrons. ux 
from Measured Activity 
The measured activity must be multiplied or di- 
Maded by several correction factorseimy order evesaceoumianod 


mOsoSeS OCCUrring In the counting precess - aceord a -are 





PAB, 


TABLE IT-1 


CONSTANTS TO BE USED IN EQUATION FOR NUMBER OF 


DISINTEGRATIONS DURING COUNTING INTERVAL 


Cross Sections* ae = 800 X 107**emé ata - 2000 xX 107“4em2 
7 
Decay Constants Ne = po ORO IOs = CrC0t97S mim = 
135952 min 
Nee oe ara 
dies unkiel 
Product of Decay A5r = 0.02489 re eh 0) 
Constants and 
Times . Not = 0.07467 al = §8,3160 
T = 5 min 
E — rs aan Actos = O21 2245 Mc nto => 13.8600 
to = Eas min 
“AgT “Ag Al 
Values of e - 0.97541 e ™ = 0.062536 
Exponentials 
“At arent 
e 2% + = 0.92805 e em'l ~ 0,000245 
eae ee 
e -° 4 = 0.88298 e OM’2 _ 9.000002 


shai ee 
e VY lle 9’ { 0.04507 


Be ual = BS = 0,000244 


*Hughes, Magurno, and Brussel, Neutron Cross Sections, 
Supplement No. 1 to BNL 325, U. S. Government Printing Office, Washing-~ 
ton 25, 0. C., January 1960. 
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+. Witten the werivity 


Cree. 1eid ) Koomuz eumere rc amd. Tay ie ic. 
beime counted is beta activity, losses aremespericnceamsauc 
to self-absorption and scattering by the foil material as 
Well-as the usual losses due po counter, €f ficiency eek ace 
peession, geOmetmyemactor, euc. 

Tf Nyas the number Of disintegrations ceceurr imme 
icine “he interval in which the néeutrone@etector 161 lars 
counvced, and Cry is the total number of “Ceunvesduring sth 


interval which are due to activity induced by thermal neu- 


Proms Only. then the re latvouship becweoens N wana Cah is 


F (ad) 


C 
N Ima Sp 


- 2.011 F F.F FL. F (a) 
yo ee) 


Ss seé&se 


fmere thesfaector Be0ll is a@ correction factor tow Gavemmaco ac-— 
eount the Maxwellian distribution of thermal neutrons as well 

as the dependence on neutron energy of the ori Viesiecesnor. the 

getector foil material. The F factors are correcti COmmeiactors 

for £amma photen counts fromesecempe time modes Cf@@ecay Inedys— 

feosium, detector geometry, counter ECLlficiéme,, Dae ae a 7 

moe. self absorption and self scatter in the fol sana 


depression by the foil. The ee and F are depend- 


sa&se 


aero DeLee Ol leathecknes's , :d.. 


*Greenfield, Koontz, Jarrett, and Taylor, IDL 
eteoutes, Jo, No. 3, ps Si eaten ie 7. 
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The totealeeoun ts Observed, CU Wid Choa emauc FLO 


rola 
aeLivity induced by thermal neutrons is Gemermamea Dye eaking 
Ome measurement, Cag, in a certain location Wiehe “heaee 
wrapped in a cadmium cover, and by taking a second measure- 
ment, C,,4,1,, in exactly the same location using the same 
tonwawi th no Cadmium cover. The difference of these meas-— 


miements is the désiredwvalue, Cy). Ihe matwend wilco <— 


Deession for Cy, is 


wien the foil material is dysprosium. For materials) viiumen 
have large resonances in their absorption cross sections for 
MeuUtrOn energies greater than thermal, the cadmium ecowewen 
foil measurement, Cg, must be multiplied by a correction 
ieee tor Bee to account for epithermal neutrons absorbed by 
moe cadmium cover... However sinsthe case of dysprosium, > 
Mmeme are no stweh large wesonances. (and thus. ssaimce C nd is 


ommsy about 1 per cent of C this correction “tsame. 


TO ve ie 


Meeessary. Preliminary data taken using dysprosium foils 


Showed that omission of this correction factor, F,4, which 


-C. W. Tittle, Nucleonics, 8, No.0 6 5 poe ines 
al. 





ey 
moa yet has not been@experimentally determinedmioer athe case 
meen the detector fo1l maverral is dysprosiun, pvoulds int, o-— 
duce an error of less than 0.07 per cent in the data taken 
mor this thesis. 

Since the beta counting apparatus used (see Chap- 
mer ILL) had 4x geometry, the correction factor for geometry 
ieee > Cgqual to unity. The scintillatver compainer dive n. 
femenris were such that all beta particles were .capturcd pre. 
eemsuncir reachitgethe walls of the contasgmer. Emie, thewecr— 


mew sOon Tactor forwback ecattering, Ff is also Cqiaite 


bs? 
ieee V 

All gamma rays emitted during the 10-minute count- 
memincvcerval are coincidental with the beta particles When 
faemerei ngs counted, so the factor ey? which compensates for 
moma Photons emitted by decay modes other than the cone being 
fered, is also equal to unity. The 15-minute waiting period 
freer tO the counting interval was purposely chosen so thay 
Pmemaunmber of gamma rays from the isomeric transition of 
ier to Dyt°= (half life of 1.25 min.) would be completely 
Mesltigsible in that less than 0.025 per cent of the ae uivitr, 


Waeesent due to this decay mode @t vime t | "O Tremaine sac — 


j= 15 minutes). 
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The= facies which compensates for the depres- 


sp? 
sion of the thermal neutron flux by the foils being used to 


Penmduct the measurements, can be calculated byethe equation® 


Pe) h yal et eee ec <i 
Mer 


witch was developed by Bothe. in this equation | K=-= sre 

radius and x = ery eT} Fee the transport, meanei wee 
1G 1+ (cos @ ave 

mom in the reactor. Ag is the Scatvering mean freee oawn 

in the reactor, and (cos ee = 2/3A where A is the mass 


number of the scattering nucleus. The term q@ depends on the 


moe thickness, d, and is determined by the equation 
a(d) = 1 -lexp (- ua)] (2 - na) + H2a?B, (-ua) 


feere G2 1s thickness of the detector toil, ff iss thewriewiae 
feweron absorption coefficient for this foil, and Es (-ud) 


is the exponential integral funetion’ value for the foil. 


Sc. W. Tittle, Nucleoniecs, 9, No. 1, p. 60, July, 
mS 1 . | 


“RE, Jahnke, and F. Emde, Tables of Functions, 
emer Pub. Co., New York, 043. 
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Calculations using this equation were compared 
Toate ex pDerimenvalk valuecemns Bothe,® and agmeacement was exce l— 
items even tor cases where the foil radius, Ryo wasmenia, 
slightly less than Ney Thus this equation was used to cal- 
culate Fsp(d) for the foils @ised Ganethis project mame 
values calculated are tabulated in Chapter LV. 

ine factor Po accounts for tbe sci ficiency Seimei c 
beta detection apparatus and is less than unity even for a 
mimecometry Getector Since the contributions by low cnemey, 
peta particles may be lost in the tube noise (see Chapter 
d) 


ie). The factor F , the correction for self absorp- 


Bane. 
mien and self scattering of beta particles by the tortoise 
dependent on foil thickness and has not been determined for 
feeeosium foils. By e@incidence measurements siti he 
fact that gamma rays are emitted simultaneously with beta 
particles in the decay of Dy?°°, the combined correction 


oc GOT Fe Featse(@)| Can Pesce eimagned,. The (equa lieme muse e 


in this determination are: 


C 


y = Kf,fpN Yor the gamma ray countves 


eel For the beta partielle ween ae 


Ce = FeFsaasclg 


af Op Tittle, Nweleontcs yes Ora ee emcee 
mS 1. 





oD 
and Ca-+ = KP oto PePsaascN for the coineidemce ecounmter where 
| Se hehe or devas ammo at Dy?°> GIsSINCesrations during wthne counting Tver — 
val and the C's are the various observed counts corrected 


meme background. Uhe factor [K f. f accounts for the geom- 


g | 
etry and efficiency of the gamma detector and for gamma rays 
absorbed prior to reaching the gamma ray detector. BE limina- 
imom Of {K fe fo] yields the equation, since Fg is kiown vo 


be unity, N = CaC,/Co_.. Aaiesere che Ce teretphe= 41 counter sis 


Pen measured, 


C_al4n counter) Co. 2 Cn 41) coun.) 
= Fe Fsa&se = eee 
N CyCe 


Since F_, is a funetion of foil thickness, d, coincidence 
se 
counting measurements were conducted (see Chapter III) for 


maem f£Oll, and the resulting values for [ie F (d)] are 


sa&se 
tabulated in Chapter IV. 

Thus, when all correction factors arewimelumoe 
Ween tinal equation for absolute thermal neutron fiux fer 2 
Biven number of counts observed during the counting interval 


ty (= 15 minutes) to t> (= 25 minutes) is given Dy euwemeqma. 


fee) 1) 


i Le eree | Sea Fey) 
th ~ {17012 Pew teawe ds) (ecm 4 2 arn 
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milere mM 1S the foil) vweleat in erams sand Pin 1S si Neu ween s 
per square centimeter per second. 
Tf the west. 18° desired. in terms on encuton 
density, n (neutrons/cm?), we must divide Dis, by the thermal 
Memiron veloci ly “om=2cu, USO centimevers per Second, Tog en 


the equation 


[Crotal ~ cal ave) 
220,000 [12011 Fe Pee. (ocean | 


were Nn is in neutrons per cubie centimeter and m i1sm@tiewo1. 


weight in grams. 
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EXPERIMENTAL APPARATUS AND PROCEDURE 


ee Reactor and External Neutron Source 


The University of Texas subcritical reactor” is 


eemposed Of a polyethylene-clad, homegeneous cOre ot 20 per 
cent enriched UOe dispersed in high strength Marlex-50 
polyethylene, and three separate reflector assemblies. 
The core is a cylinder approximately 10 inches in diameter 
ond approximately 14 inehes lone, assembled tromeerten 
homogeneous polyethylene discs impregnated with UOs. One 
radial and one axial hole of approximately 1 inch diameter 
memebrate the core for the purpose Of inserUtenso, owe, 
MemLrOn sOUurces, etc. 

The reflector assemblies are made of polyethylene 
and graphite. The polyethylene reflector is composed of 
iorsconcentric shells, each approximately three — mene. 


tadck, sO that a total reflector thickness Of €lcie mame: 


°F. D. Eichenbaum, Lockheed Nuclear Products 
Reports NR S?7 and 98, March, 1959, Nuclear Prodwe sare. 
eon Of Lockheed Aircraft Corp. , Martetta, Georeuae 


on 
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or six inches may be used. The graphite reflector is ap- 
proximately ten inches thick. All three assemblies are 
Constructed in the f£6rm of annular cylinders with end plugs 
to fill the annulus so that the core can be completely re- 
flected. 

During experimentarvion, the core 16 inset ceee co 
one of the reflector shells and the appropriate plug is 
meserved into the end gap, thus providing a uniform vier 
Mess Of reflector at the ends and around the core. When 
mime cote is properly positioned, its axial and radial peace 
are aligned with similar ports provided in the reflector. 
ims Tolls, control rods, and neutron sources; apprope sac 
tes vie experiment being conducte@, can be pesittiemed in 
tae ports. 

pmall.s cylindrical plugs, scme compoced eo! care 
feeserial and seme of reflector material, eare=provrdedanin 
aeerecient quantity fOr filling, as appropriate  toeunc 
epeepertement beingeconductéedyethe yvotds createdein the re- 
actor by the radial and axial access ports. 

In this experiment the 6 inch polyethylene re- 
flector was used. A cross sectional view of the core with 
the 6 inch polyethylene reflector in place is shown in 


Higure ILI-l. The dimensions in this» figure aremm 





oe, 


Gentimeters and are accurate dimensions rather than the 
approximate figures used in the above description. Dimen- 
exons of the individual core element discs as well as the 
various constants applicable to the core and the reflector 
are tabulated in Appendix A. 

The cadmium control rod used in this experiment 
mes & hollow cylindrical tube 86.a centimeters longeicog@- 
posed of commercially pure cadmium. This tube had an out- 
side diameter of 2.62 centimeters and a tube wall thickness 
of 0.09 centimeters. 

Pres boron coOnurol red. was sane lc yee nde | 
tube 80.4 centimeters long composed of polyethylene 4 per 
cent impregnated with boron. The outside diameter of the 
tube was 2.71 centimeters and the tube wall thickness was 
Ojl5 centimeters. 

Because the University of Texas Subcritical Re- 
meemer 1s subcritical, an external souree Of Meuvremo musL 
ime used in order to maintain a thermal neutron flue gwieain 
the reactor under steady state conditions. Several manu- 
factured Pu-Be sources with a constant neutron yield are 
available for use with the reactor. Unfortunately Sade y-— 
ever, the dimensions of these "fixed" sources prohibit 


their use inside the reactor when one of the control rods 
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PoopOsttLOned in the axial hole Thus. uy Was neceocary 
to use the Van de Graaff accelerator to create the neces- 
sary external neutron source. 

The Van de Graaff particle accelerator available 
for use in connection with this experiment was manufactured 
by High Voltage Engineering Corporation and is capable of 
producing deuterons accelerated to energies of up to l Mev. 
A beryllium target was bombarded by l Mev deuterons from 
the accelerator and neutrons were produced by the 
Eee (din) B+? reaction. The Location of this external neu- 
Gron source within the reactor is shown in Figure LIIi-l. 

Beryllium was chosen as the target material be- 
Simseesuhe sSsuabilitcy characteristics Of the accelerator 
required operation at deuteron energies between 0.8 and 
1.0 Mev, and, in this deuteron energy range, the Be?(d,n)Bt° 
reaction afforded a better neutron yield than does any 
other X(d,n)Y reaction. Also, since beryllium targets 
were avallable inmemetal discw@ioerm, the probilems aseeciacved 
with target cooling and handling were greatly simplified 
by use Of beryllium rather than tritium, the oOnlyaawucn 
target material immediately available for use in connec- 
tion with this experiment. Target cooling was accomplished 
by means Of a stream of air blown directly on the target 


assembly. 
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The energy distribution of neutrons from the 
Be?(d,n)B?° reaction for 1 Mev deuterons (as shown in 
Figure III-2) shows a preponderance of neutrons having 
esmergies Of approximately 3.9 Mev and 4.7 Mev.” “Phoevencu- 
lar distributiont®? curves for the neutrons produced by 
mimes reacuton indicate only slight directional pretereneccy 
and measurements conducted to determine the strength of 
mie sSOUTrce Used in this experiment showed Chav. if san 5ecr 
fropi1c neutron yield were assumed, the maximum deviaron 
due to directional preference would be less than 6.5 per 
gent. Thus, if Was assumed fOr the purpose Of Lites semper 
ment that the Be?(d,n)B*° reaction functioned as an iso- 
tropic source of neutrons. 

For this experiment, the source strens ieee 
be determined as a function of the electrical current cre- 
ated by the action of the beam of accelerated deuterons 
Peoduced by the accelerator and striking the beryllium 
target. This was accomplished by measuring with a "long 
counter’ the number of neutrons produced at various beam 
current settings. These measurements were taken at various 


angles relative to the deuteron bean. 


1°A. I. Shpetnyi, Soviet Physi cometh omens 
Mewober, Llgs/. 
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N(E,) (Summed over 
all angles) 
ie = ne 8, Mev 
E, (Mev) 
O al la 2 4 SD 
mreure IIll-2. Energy Distribution wot Neu Cron cae eon 


Be?(d,n)B?° Reaction.--From Shpetnyi.* 


*A. I. Shpetnyi, Sovaet Physics JETP > Oe ee 
Bemober, 1957. 
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A calibrated Pu-Be isotropic neutron source with 
a yield of 3.5 x 10° neutrons per second was then placed 
in the same location that the beryllium target assembly 
previously occupied and the neutrons from this "known" 
source were counted for the same time interval using the 
long counter. The average of the measurements taken at 
various angles using the Be*(d,n)B?° source was assumed 
to be the measurement due to an “equivalent isotropic" 
mMeutron source. This average value was then divide@eby 
the measurement taken using the calibrated source and the 
resultant quotient was multiplied by the strength of the 
Calibrated source in order to determine the Be” (dja 
source strength in neutrons per second. Thus, the "“equiva- 
lent isotropic” source strength for various deuteron beam 
currents was determined. The results are shown graphically 
mie tgure JTT-3. The source strength for a beam *currene 


of 1.5 microamperes was found to be 92.2 X 10° neutrons 


mer second. 


Be Foils 


As previously mentioned in Chaprverel! Seen 
herent operating characteristics of the Van de Graaff 


accelerator used to create the external neutron source 
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dictated the use of dysprosium for the neutron detection 
Magemial. The stable: Dy--- icotope weetul no tht eee i 
ment comprises 28.2 per cent of natural dysprosium. The 
Other isotopes contribute essentially nothing to the activ- 
iy measured aGuring this project. 

Since no information concerning the self-absorption 
Characteristics of dysprosium could be foOUnR@@in the Titgera— 
tuie, the variation Of beta counting rate Wich) for cher — 
ness was determined experimentally following standard 


procedures.t? 


Dysprosium foils Of several different thick- 
Meocoes Were irradiated in the subcritical reactor at tne 
geme location for identical time intervals. After irradi-— 
en Of each foil, 15 minutes were allowed to elapse pmuer 


165M yould vanish. 


to counting so that the activity due to Dy 
Bach foil was then counted for a period of LO minutes. The 
results are shown graphically in Figure III-4. This self- 
absorption curve shows a@ broad maximum for foil thickness 
between 90 mg/em= and 130 mg/cm. This corresponds to 

Moe kness between ©.0055 inches and 0.0065 inemes, “comune 


foils were fabricated with thicknesses as close as possible 


momo tLnousandths of an inch. 


tig. W. Tittle, Nucleonics, 8, No. 6, P- 5, June, 
BoD] . 
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Thousands of 
Counts per 
LO Minutes 





0 ao) 40 60 SO 100 i120 140 eo Petree 


Foio Thickness (mg/cm*) 


Pore TTI-4. Variation of Beta Counting Rate with Moar 
Thickness for Dysprosium Foils pig oicmeernd 
Diane wer. 
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The foils used in this experiment were fabricated 
from a sheet of commercially pure natural dysprosium metal. 
Pieces of the metal were turned between centers on a lathe 
to a diameter of 1 centimeter and were then ground by hand 
wor a thickness of about 6 theusandths of an inch. 

Upon completion of the cutuing and grinding a 
foils were measured as accurately as possible. The foil 
weight was measured on a Wm. Ainsworth and Sons arm balance 
to an accuracy of 0.1 milligram. The diameters were deter- 
mined by averaging five measurements taken with a micro- 
scope designed for measuring tracks made by particles in 
mmeLear emulsions. The area was then cemputed, ancmine 
foil thickness was determined in units of milligrams per 
square centimeter. The results of these measurements and 


Sarculations are tabulated in Chepiter viv. 


to Beta Counting Equipment 


® 


A liquid scintillator, composed of 4 grams per 
liter p-terphenyl and 0.1 gram per liter "Popop" dissolved 
moecOluene, was used in Order to Create 4x geometry ane. 


beta detection.+= This scintillator was contained in a 


t-p, T. Doss, Unpublichedethesis, Uniyer- so. emom 
ewes i Ute jeelg6 4 
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eylindrical aluminum sleeve with a removable aluminum cover. 
This sleeve (see Figure III-5) screwed into a machined 
brass sleeve which covered the 6810-A photomultiplier tube 
used to detect the scintillations created by passage of 
beta particles through the liquid scintillator. 

Two teflon "0" rings were used to create a seal 
between the photomultiplier tube and the aluminum sleeve, 
mpiues holding the Liquid scintillator in the upper chenbere 
The removable cover, which provided easy access to the 
liquid scintillator, had to be designed so that the system 
was entirely light tight when closed. This was accomplished 
by means of a concentric series of "tongues and grooves" 
machined in the surfaces of the aluminum sleeve and the 
cover. The foils were suspended in the liquid scintillator 
by means of a small spring clip attached to the removable 
cover. 

The dimensions of the aluminum sleeve were such 
fae the amount of liquid scintillator surrounding sinew ior! 
was sufficient to absorb all the beta particles emitted by 
Bao decay. The maximum range of the most energetic 
(1.3 Mev) beta particle emitted by Dyt©®> was determined 


by the egquationt®? 


13s, Glasstone, Sourcebook on Atomic Energy (2d 


ed., FPrineeton, N.JdJ.: . D2 Vanier ond, ver, Lose Mee lew 





Aluminum 
Guide Rod 


et von 
O Ring Seal 


Proere Til. 
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st a 


Aluminum Cover 


Aluminum 

SC i Gielelest our 
Coptainver 
Sleeve 


Proud scinti 1 avo 


© 
6OHOre Phieto— 
Nutctiplier supe 


Brass Sleeve 


Details @f 45 Gecmei ry sbi¢dui deci 
tillator 86a Deteoevs. 
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Ei. ~ O.17 
Meee (2) em) = See es BG. O.5 Mev < Emax <c 5 Mev 


lees 

Or Bay = L.S Mev, this equation gaye a Maxinumerance or 
0.595 g/cm=. Toluene has a density of 0.866 g/cm®, so 
Rmax Was 0.688 cm in the liquid scintillator. The smallest 
Batekness Of scintillator adjacent to the foil was 0.750 
eae. Gli beta particles were absorbed prictm.e Peach ene 
the walls of the container. 

The 681OA photomultiplier tube was mounted in a 
Mexas Nuclear Corporation tube base utilizing a 417A cath=— 
Ode follower. A block diagram of the entire beta detect- 
ing and counting apparatus used in this experiment is shown 
in Figure III-6. 

The proper discriminator setting, fTometne ybaed- 
Atomic amplifier was determined by increasing the setting 
from zero to the point where the observed counting rate 
approached the background level. At this discriminator 
setting, the photomultiplier tube electronic "noise" was 
Serminated. Undoubtedly, detection Of many Gi ine ven, 
low energy beta particles was also eliminated. However, 
mois loss is corrected by the FuFoagccl(d) factor previously 


discussed in Chapter II. 





DC 


I9TeOS 


o&l LIPONW 
vTWoLTyY pared 


ol i 


4096 TePOW 
OTWOLY pateg 








Beigel 


-eieddy Sutzunoy pue uotyzo9490q S4aq JO wWerseTqg Yoo tT, "O-TITI sinsty 


ATdane 
| IoOYVeTnNsesey LIMOd AH 
ie SET VY G2T2 latofo 09 Sly Tepon 

Tepoym ussuasog Iauuey 





SSCL TOA USTH 





SToOAD Q9 
SICYTOA MOT 
| ISOMOTTOT 
See LaTITTduy spouyeo VLTV 
7np SUT PBOTIZAOUON] | ULITM 





O71 B68nNUs _ 
' maar! GTZ L9PpoW 


oTWOLy-pated 


\ 





sseq sqny, °dao0g 
IBSTONN Sexo, 


94 NY, 
AST OC a ey 
-070Ud 
VOTS9 


ROLO TT Tmt | 
Din Ome 








oS 


The factor F.Foaggc(d) was determined experimen- 
tally by coincidence counting techniques. The various 
foils were irradiated long enough to reach saturation f6r 
tHe 2.52 hour half-life activity wt Dy*°-. Fifteen @emmercs 
were allowed to elapse between the time of removal of a 
fiom from the irradiating flux and the™ time of counvine. 
Weroil was placed in the beta counting apparatus and then 
a gamma ray counting apparatus utilizing a NaI crystal 
gamma detector was positioned adjacent to the aluminum 
sleeve which covered the plastic scintillator for the 
Reta counter. Gamma rays emitted coincidentally with the 


Bema decay of Dy*°? 


passed through the aluminum sleeve 
Smamwscome were detected by the Nai-crystal. 

me cOlmemagence CoOmmir, CR ys was measured first. 
ire, gamma ray count, CA, and the beta particle counts, CB 
and Cp (4a counter) were then measured for exactly the 
Same time interval, and the observed Cy and Ce va Lue 
were projected back along the proper decay curve on semi- 
ieee sraph paper in order to»normalize the Cy, Ce; and OBL 
counts to the same time interval. The factor F.F,,¢cce\@) 
was then determined by dividing Ce, x Cal4n counter) by 
Cy x Ce as discussed in Chapter If. The values Of 
ihe 


= areeilids) thus determined are tabulated with other foil 


taco Chepiter say. | 
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D. Experimental Procedure 
The foils were exposed, one at a tine. at yvariees 
locations in the radial access port of the reactor assembly 
at intervals of approximately 1.d@ centimeters. The irradi- 
Qimon tweme for each «fedl was S minutes. After waltimemisd 
Mm tbes after removal of the foil from ™cne reactor, the 


165 was counted for a 


bewa activity from the decay of Dy 
period of 1O minutes. During the 15 minute waiting period 
a 10 minute background count was taken and this value was 
assumed to be the background count during the foil count- 
ing period. 

This procedure was followed for both the cadmium 
covered foil measurements and the bare foil measurements. 
The same foil was used for both measurements at a given 
iereatvion. During the bare"-forP measurement, a pelvetnyvene 
Spacer disc was used to replace the lower cadmium cover in 
Order to position the foil in exactly the same location 
iiserar as was feasible. The actual difference in location 
was never more than 0.53 millimeters, and since C,,g is only 
aim@ut 1 per cent of Ci;, the error due te any difference 
in location was assumed to be completely negligible. The 
cadmium cover used was in the form of a cup with a disc 
wadeh fitted completely inside the lip. The cadmium used 


moYthis cover was 0.051 centameters tuick. 
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The beam current (and hence the source strength) 
for each irradiation was determined by taking 2O beam cur- 
Memy readings, equally "spaced tn time, duving hewn aioe 
irradiation period. The average of these 20 readings was 
wie value used for the beam current. Al] Poil activity 
measurements were normalized to the same beam current of 
imo mMicroamperes. This amounted to normalization te an 
external source strength of 92.2 X 10° neutrons PeResecoud., 
since the source strength is directly proportional to the 
beam current, and since 92.2 x 10° neutrons per second is 
the source strength corresponding to 1.5 microamperes of 
deuteron beam current (see Figure III-3). 

The scintillation counter was kept completely 
fol of liquid scintillator at all times, and theyligquad 
scintillator was changed frequently in order to prevent 
Peiicaminatlion.s ihe folls were Kept as clean ac. perc tale. 
and were washed in used scintillator each time prior to 
being inserted into the detector. 

Prior to commencing measurements each day the 
counting system was checked for proper operation using a 
Na“= source of constant activity. This source gave a 
@ounting rate of c8090 counts per 2 minutes. Thesvacite. 


equipment settings used were as follows: 
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High Voltage eee LIS IOLO Migoml ine 
Ap Waiteer Ubsoecriminator =. Vee. poo. 
Amplifier Cain» « «-. « ».@ 76260 
Scalar Discriminator 17... “On 


Measurements were taken first with no control rod 
Mmieune reactor and Ghen with the cadmium control rodmainerne 
ieachor, and lastiy with the boron control rod in ineuce 
meEpor. The control rods were placed along the longi tugqgiuna 
axis and covered the entire length of the axial access port 
Of the reactor. The effect of the boron rod was semen 


Was not reported in this thesis. 
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RESULTS 


The radial. Neutron density 7dis tert puro «cul ees 
pee LOrmined "by "experimen talkie ra cds vad ome Om mine 
wee IV-1l for the reactor with and without the control rod 
ie position along the longitudinal axis of the assembly 
Mme neutron density values showm here nave been “emma med 
imo constant source Stremeth of 10° neutrons per second 
from the Be(d,n) reaction. Comparison of the two curves 
auewe oa large depression of the meutron dems hiss clic 
feeinity of the control rod and a general decrease throughn— 
Gime. the reactor. This 2s to) be expected s tteee ie aeaigc 
memoer of neutrons absorbed by the cadmium rod are no 
Meeser available to diiiuse through the reactor Of Steere 
immeipate in fission events whieh release other neutrons 
wemcvhne system. 

The single veloe@@mty groupedaifusion theory sr een. 
mmpue was used to prediet Chew mweutron density, dis Uo m 
for each case mathematically. The results of this treat- 
tem. are shown graphically .in  eweure [V—-c for thesreseran 


with no control rod and in Figure IV-3 for the reactor with 
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me control rod. The agpeement between themthecore tiga 
and experimental distributions could be termed no better 
ima Tair in E€ither case. 

The large "hump" observed experimentally at the 
boundary between the core and the reflector is not ade- 
meavely represented by the theeretlicalmeiet, soutien. al 
though a slight change in slope is indicated at that posi- 
moon. it is interesting tomtiote, however a"ha hee orca 
under the theoretical curves are essentially the same as 
ihe areas under the corresponding experimental “eunves. 
Hous implies that the theory does provide “Torewtrespnepes 
botal number of neutrons althouen the distriputrrone. 
somewhat distorted. 

It is believed that the agreement between the 
distributions predicted by diffusion theory and those ob- 
served experimentally would have been better if a more 
meeurate value had been available for the bucklingewoit Suaec 
meme tOr used in this experiment. Calculations shovegnuat 
if the value Be = 0.02cm™~* were used in this case, the 
MmamorOn density Predicted einconepuca ly LOR unc secraan 
reflector interface would have been almost exactly the 
Same as that experimentally observed. 

Analysis of the davavwobptained for varioueepe «2 


tions at the core-reflector interface indicate that surfaces 
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of equal neutron density would have the form of an ellip- 
soid of revolution with a ratio of major and minor semiaxes 
Pmopeproximatvely, 1 csvom i. ris (mp irvesr that . ie ine 

eoamee in slope ab thevcere-retlectvor intertace Verecemece 
ieeved, theoretical calculations for the system mignursoe 
made with considerable success using spherical geometry. 
Bomculations involving spherical coordinates are Sommcemm.—. 
Heh more simple than are these involving cylindricaweee. 


Seadinates. 





AVP PE Ni eiaie ox A 


REACTOR CONSTANTS AND CALCULATIONS 





bo 


Core Dimensions* 


Fuel Hole 

Element Thickness Diameter Diameter Volume 

Number (inches ) inches inches cubic inches 
al 12900 SSMS) Le2ee 144.18 
Z Ze0D0 Si. S10re! 12286 Le Seae 
3 ya Oe, 927306 1 200 156.84 
4. 2207-0 9.905 1260 156.84 
S 2.045 9.912 1.260 WSS in? 
6 23051 9.934 1. Ze2 LSG276 
7 ieOoo oe 2oL i 2318 78.84 
8 NOT USED 
g OF 507 92 2Ze 1.285 SGaoe 

LO NOT USED 


Total Volume = 1,042.54 


Expressed in Terms of Cubic Centimeters, Centimeters, and Grams. 
Total Volume = 17,090.36 cm® 
nverage Diamerer = 92.914 inveae2s-1oc em 
Wotal Thiekness = 15,754 ini] 54.66o cm 
Total Weight of Core = 18,603.22 grams 
Core Density = 1.08852 grams /em® 
Total Weight of U0? in Core = 2,668.49 grams 
Total Weight of U°°° in Core = 465.04 grams 
Total Weighv or Uo dim Core = 2,55) 054 grams 


*Data furnished by manufacturer, Lockheed Aircraft Corp., 
Marietta, Ga. 





70 





‘gcet ‘Aine 6°09 ‘a ‘Gg uoqSuTysem ‘a0TTIO SuTqutag “yon ‘S’n “See "INA SUOCTJOSeg SsozQ UOTINEN 
‘-n oy Szqremuog pue *‘°p °q fsausny Worzy uayeq, are pue suorT{NSU [eUIsy, TOS aie suOTJOSES SSOLD TIV x 





9070°O 26le°T © (ode © 2520 0 Naive o/TV.LOL 
Z2TTO°O 0000°O O£00°O 8cg00°O Of00°O 8s000°O0 2200°O S000°0 eG 28 og0l * V6L°S — geal 
TéT0°O 9070°O *7870°O 8200°O 78700 L270 0 2000 0 0° 28S O°o89 O°OT oe0T * 869°0~ egal 
GO?ST*O O 8200°O LT¥0°O 6200°0O 0000°0 6600°0 0 2000° 0 o°V = =4e0T ¥ P86°9 9t O 
2LST°O OF SO¢el 0 es 000 L26T°O TOOO"O 9dé6T*O 0 %E00* 0 8° s50T * G2eTO’? atv 
0000 °T O  O870 1 #999 6 ~ 1970" oe 20ce Osco 0 eae, © O°8S 3s50T * Gd0'8 iH 


a A 


1.0 e S e S 
a ak BAX 9s00 ) aK EG < FO 0 0 N 


eh pr CP A 


(wo/us goeso’l = Aytsuep) xsqueqsuoQ e109 





(ele 


BSio nie 2a, ak 
‘og ‘uoqgSutysem ‘aoTsJO BuTquTag *4A0D *S'n *GZEe INA suOTJOSeg ssorM UOTYNEN *°T *Y 
‘2n.remyog pue *‘*p *q ‘Ssausny Worly usyey ale pue suo NSU TeULTISsyy AO ate suotzossS SSOLO TTV x 











Slo’ 1 CVOL’ ZS 999¢°S G/Z0°O T8Ee'"s STV.LO“ 
VLG Orie oO 9SS0°O T86L°O TOOO’°O E¢86T’O 7VE00°O 8° 7 avOl * SITELl’? atv 
6EL10°T Cob60'? 1999°0 CLOL’S PL20°O BEETS oos O O’°es avOI X 9292°8 iH 
VE Sz_4dee wo ro /9 = wo wo wo surec SuUIeC wo /su01e 
(Eee) (im) (we/e=) (pm) (G8) (qm) (SuHeR) (sured) — (guo/smoae) 
¥o/2 92 =" "(oc0>) eke < ax eK Po =0 N 
(_-wo)t4¢ 


—— 








( uo /us L66S6°0O = Aqtsuaep) xsqueysuog 1oJdaeT ey 





CZ 


REACTOR CALCULATIONS1# 


The transport Mean iree pathy), 9 ens auetic ier ies 


Given by the equation 


ile IL 


X —— 
ieee = ( ] ) 
Ltr atotal Lig 2/ SA 
where 
bt 


: E 2a (C2 Bee) 
r = *total 2 Cn -» eee are ZA. 
il 


fee one i elements of Which the medium 1s compoced tau. 


Memeeuie reactor core, 


(c) i il 
ee 1.2792 em + 


By the same method for the reflector is 
af 


(R) it iL 
Xr = ee eS ie : 
aoe “TR) O 792835 cm 
z 1.2613 em 2 
J 


14Meghreglian and Holmes, Reactor Analysis (New 
monk: McGraw-Hill Book Comgmoaojpmcnep. 05. 





i 


ive diimusiom coethicient, Do. tor aemed iin creased 


im; the equation 





Piere At, 1S the transport mean free pevun. Om) Pile eceree D. 


is 


D = ee = Ov ok 74~em = 0.26058 em 


mie cor the refleetor , D. is 


D = teat eS Oo coon em = 0.26428 cm 


The extrapolated radius p, at which ¢ becomes zero 


Ls 


P59 = Pmeasured + 2Dp = 28.32 cm. 


The extrapolated height z at which ¢ becomes zero 


is 


Zo = 4measured +2 DRO =" 55.55 cm. 





74 


The values of %p end 2 at sone 1ntvertace spe Weenie 


core and the reflector are 


1S oS Saien 


@ 
ye 
I 


and vd ly 244 een. 


The diffusiom length. Eo isedefinedes, Sele sega. 


Miewet- = D/x. Thus, for the come, 


p2 - DB. _ 0.26058 em _ 5 yee one 


¢ x(c) OU SE Seimee 
a 


and for the reflector 


ee = = Os eee meus Se Oenl@ em< 
5 00275 tena 


a 


The constant KP is defined) by. the. cq ua tio nee 


2 
R 
ew Thus; Lor the (remiecier, 


In Chapter I it was determined that ie n/2z0. Piet bisy ye aian ie 


the reactor used in this experiment 





(e 





I —i 
a — = 0.0468 cm 
Z a 
Since «= is defined to be n= = k* - Ro tier Germ 
R R O Z 
Ky can be determined as follows: 
Ky = (eo 2 aoe CRONOE mn OnOOcIo) — 0. liaie cm. 


In OF¢@er.loeC¢ cecum nice. vie wequa Tron 


Gain ees Z = z ane E Z 
Za eiZie, Bh CZ 6 ue 








Ke 
tire be solved by numerical methods. This solution syietd. 


a value of 


ky = OG Apecim = 


Using the value wou Ba furnished by the manufacturer 


@eeone reactor, k can now bemGetermined by thewequs tied 


ors 


2 2 
Ka) Beis 
Since B= = 0.027172 em~~, this relationship gives the value 


of 





at oe and 


78.7 n/em® 


Powlows: 


giving CR 


cS 


Ky = 0.1580 cm 1. 


Evaluation of ¢COnstants involving Bessel Funeticns 


by are as follows: 


Kono.) EKGs 


TC = =—en Oo 25 
LAR a (GRASS | 
O~ SOE O 


b. = b = 6€ =] 902795 3c 6a rea eee 


Jo (kate) Se Oae are) 


a 0.465. 
UE De) ¥ (0.0579) 


For the case with nO "con vrel Ted ew ae Jae = 


when p = z= 0, the constant Cp, is evaluated as 
Jylgke = 
Cr = Ac o (kp; ) 
E ( ) Ko (pp) eet } 
een) (ele 
O Oma Ig (k,00) OMe ec wail 
J ee ss) 9 
eS 7 gi n/cm® 


[K (1.4371) - 0.0045 I9(1.4371)] 


= 80.56 n/em®. 


When the control rod is in position, A, is taken 


to be the maximum experimental value in order to correlate 


experiment and theory. This gives. the wal ue@iee—=>.. 2 n/em>. 





fig! 


Berne this value, Cp can be evaluated by the relationship 


J_(k b_) 


O O 
aS 5 Cp0.) " fo (rg01)] 
ie Cc K Ck } 
7 ome . 
Fo egos) Tt (r Fols304)] 


Po) 


[J (1.989) + 0.465 Y (1-989) ] 


54.3. — ____ fem 
[K, (1.437) - 0.0045 I, (1.437) ] 


Biving Cp = 1llz2.73 n/em*. 
In summary, the equations for n(p,z) with all the 


constants evaluated are as follows: 
(a) Without control rod: 
i), 2 ) = 78.7 J (0.158) Coss. 04/14 
Deon 2) = 80.56[K, (0.1142 ) - 0.00451, (0.11429) ] 


x COs Oo, U4Gez 


(b) With control rod: 


n,(o,z) = 54.5[J, (0.158) + 0.465Y (0.158) Jeos O47 


np(o,z) 112.73{K, (0.1142) - 0.00451, (0.1142p) ] 


“% cos (02046672 


tmese equations were used to determine the theoret wea enews 


density distribution curves shown in Figures 1)-2 ana —o. 





A PF PEN D 2x 


PROBABLE ERROR 


B 





is, 
The probable error of a function F(x,y,z) is given 


Ey 


Pe. eg OE cin 
OX Y oy 4 Oz 
Wier e Ps Be and P. are the probable exrors ine, eamceez. 


meopectively.  Thesprobabie eirons PL» oe and re Bae Ten ned 


fae tne general equation 


td 
| 


= 0367 250 


where o is the standard deviation of the value observed. 
The function, n(p,z), with which this experiment 
meamconcerned, is revated to thevmeasuvede quantities. Dy suc 
equation 
KC heeds 


ni ogee = ___)) SS) 


ma Fe F atac \?) 


where K is a constant having the value 


1°H, Margenau and G. M. Murphy, The Mathematics 


Seeebysics and Chemistry (2d ed.; Princeton, N.J.: D. Van 
Neoetmand Co., 1956), Sections “Waes2—o56., 
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: og 
= BE xe 


3.2 x 10° <* Seo412)< Peo cee ao 


when source strength of 1 xX 10° neutrons per second is used. 


The probable error for n(p,z) can thus be written 


Boe (1): P. = K 2/P ae 
z th nF, sare Cio et aa ae 


F wale 
an ep Ss ” an sp 
oe Fsaase eT ae 
Foe Fgaase m Fe Fgaasc 


assuming K has no error. 


= 





The data obtained experimentally for thes point 
oe 12.6c CM, Zoos 00 will Demwsed Tomdetcrmimemiie. pine bane 
erzor for this experiment. 


fhe value Cl, 15 jobeained = by Stic Trelaeaen 


Cth = CLotal - Cog 


where Ci,1,] and C.qg are measured quantities. The rea ble 
eiagor P, Wes) 
fe al 
P = ¥ P& + PA 
Cth CLotal Cod 





oul 


where 
pe = (0.6745 V17?475)|- = 6iibed 
Ctotal 
and 
p= Or 6745 402)- —a6400 
C 
ed 
Thus 
P = \/81134 + 640 = 286. 
C 
th 


The probable error in Pep (a) a8) alee ay and 
Sp 
mee Of the average error of 1.5 per cent tabulaved boy 


Hictle+° for the function F (ad) gives Pr =| ©. oles S. 
sp Sp 


The probable errom inwthe foal weigh? eee 
Omoo001 grams. 
The proba bile ~ cine eye yn for the self absorp- 


e sad&se 
tion and self scattering correction factor is determined 


Dyeune relation 


+°c. W. Tittle, Physical Review, 80, pl i7ecummooon 





She 


P = £ 2(P + (P Bay 
For saase Cay Ce C4 Ca (4x) a) 


“p-7"p'* Spay gl4™) \ 
F (=, Ee ga Fo, 2 
C Oe 


since F_F = { Poe mere f 15) aeconevent having 
e sa&sc Cec, 


mae value £ = 0.896¢. Using Veuw@e valves 


ails = 0.6745 (Cy 
im the above relation gives 


Po - =) O,00CeC66, 
e sea&se 


Thus, the total probable error 13 dere rmiiecmce 


be, using equation (1), 
P, (experimental) = 0.466 n/em® 


mest Of which can be attributed te the Jew Cg-y count rates 


(d) 


experienced when experimental determination of F F 
e sak&se 


wes accomplished, and to the prepaptewer aera Fay: 
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This value represents a total probable error of 


P, (experimental) = 1.77 per cent 


for the neutron density valwtes determined in this expen ene meee 
The probable errors in the constant Verma sus dea 

the development of the expression for neutron density are 

those give. by the authors making the determinations and 


aac. 


P Aone LOO barns tom co SOO Darnme 
o ny 
ny 
4— 
P 4-5m = 200 barns j o7, ony =o) ClO maieait Wice 


ny 
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